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Numerical simulation studies in 2D with the addition of noise are reported for the convection of a 
supercritical fluid, 3 He, in a Rayleigh-Benard cell where the fluid parameters and cell height L are the 
same as in published laboratory experiments. The noise addition is to accelerate the instability onset 
after starting the heat flow across the fluid, so as to bring simulations into better agreement with 
experimental observations. Homogeneous temperature noise and spatial lateral periodic temperature 
variations in the top plate were programmed into the simulations. A speed-up in the instability onset 
was obtained, which was most effective through the spatial temperature variations with a period 
of 2L, close to the wavelength of a pair of convections rolls. For a small amplitude of 0.5 /iK, this 
perturbation gave a semiquantitative agreement with experimental observations. Results for various 
noise amplitudes are presented and discussed in relation to predictions by El Khouri and Carles. 

PACS numbers: 44.25. +f, 47.27.Te, 64.70.Fx 



INTRODUCTION 



In recent papers, convection experiments of su- 
percritical 3 He in a Rayleigh-Benard cell with a 
constant heat current q were reportedQ, Af- 
ter q is started, the temperature drop AT(t) across 
this highly compressible fluid layer increases from 
zero, an evolution accelerated by the "Piston Ef- 
feet" 0, B B. Assuming that ^ laIg e r than a 
critical heat flux necessary to produce fluid instabil- 
ity, AT(t) passes over a maximum at the time t — t. p , 
which indicates that the fluid is convecting and that 
plumes have reached the top plate. Then truncated 
or damped oscillations, the latter with a period t osc , 
are observed under certain conditions before steady- 
state conditions for convection are reached, as de- 
scribed in refs.P,Q]. The scenario of the damped 
oscillations, and the role of the "piston effect" has 
been described in detail in refs.jg and |t| and will 
not be repeated here. The height of the layer in 
the RB cell was L = 0.106 cm and the aspect ratio 
r=57. The 3 He convection experiments along the 
critical isochore extended over a range of reduced 



temperatures between 5x10 4 < e < 0.2, where 
e = (T — T c )/T c with T c = 3.318 K, the critical tem- 
perature. The truncated - or damped oscillations 
were observed for e > 0.009 and over this range the 
fluid compressibility varies by a factor of about 30. 

The scaled representation of the characteristic 
times t osc and t p versus the Rayleigh number, and 
the comparison with the results from simulations has 
been described in ref . |8j . Good agreement for the pe- 
riod ioscwas reported. However a systematic discrep- 
ancy for the times t p shows that in the simulations 
the development of convection is retarded compared 
to the experiments. This effect increases with de- 
creasing values of [Ra corr — Ra c ] , where Ra COTT is the 
Rayleigh number corrected for the adiabatic temper- 
ature gradient as defined in refs.Q,!^] and Ra c is the 
critical Rayleigh number for the experimental con- 
ditions, 1708. This is shown in Fig.l of ref.Q, in 
particular in Fig. lb) for e = 0.2 and q — 2.16xl0 -7 
W/cm 2 ([Ra c °" — Ra c ] = 635), where an experimen- 
tal run is compared with simulations for the same 
parameters. Here clearly the profile AT(t) from the 
simulations shows the smooth rise until the steady- 



2 



state value, AT — qL/X — 95 /iK has been reached, 
where A is the thermal conductivity. Only at t w 90 
s. does convection develop, as shown by a sudden 
decrease of AT(t). By contrast, the experimental 
profile shows a much earlier development of convec- 
tion. Fig.l of ref. 8] is representative for the observa- 
tions at low values of [Ra cmi — Ra c ]. At high values, 
both experiment and simulations show the convec- 
tion development to take place at comparable times, 
as indicated in Fig. 5b) of ref.[8j, and specifically in 
Fig.2 a) of ref. |a, where [Ra co " - Ra c ] =4.1xl0 5 . It 
is the purpose of this report to investigate the origin 
of this discrepancy by further simulation studies. 

CONVECTION ONSET CALCULATIONS, 
SIMULATIONS AND COMPARISON WITH 
EXPERIMENTS 

El Khouri and CarlesQ studied theoretically the 
stability limit of a supercritical fluid in a RB cell, 
when subjected to a heat current q started at the 
time t — 0. Their fluid was also 3 He at the crit- 
ical density, and the same parameters as in ref.Q] 
were used. They calculated the time £i ns tab and also 
the corresponding AT(ij nst ab) for the onset of fluid 
instability and they determined the modes and the 
wave vectors of the perturbations for different sce- 
narios of q and e. For t > imstab inhomogeneities in 
the RB cell and noise within the fluid will produce 
perturbations which will grow, from which the con- 
vection will develop. An indication of the growth 
of convection is a deviation of the AT (4) profile in 
the experiments or in the simulations from the cal- 
culated curve for the stable fluid (see for instance 
Eq.3.3 of ref 6]). It is readily seen from simulation 
profiles such as Fig. la) and b) in ref. 6] that the de- 
viation becomes significant for t only slightly below 
t p - the maximum of AT(t). In simulations, the 
effective start of convection can also be seen from 
snapshots in 2D of the fluid temperature contour 
lines at various times, as shown in Fig. 5 of ref. jiol] . 

P. Carles 11] has argued that the reason for the 



discrepancy for the time t p between experiment and 
simulation is that in the former, the physical sys- 
tem has noise and inhomogeneities which cause the 
perturbations beyond iinstab to grow into the devel- 
oped convection. By contrast simulations have a 
much smaller noise. Therefore in the simulations 
the perturbations take a longer time to grow than in 
the physical system, leading to a larger t p than ob- 
served. Carles' comment led us to try as a first step 
imposing a thermal random noise on the top plate 
of the RB cell, which was to simulate fluctuations in 
the upper plate temperature control of the labora- 
tory experiment. The temperature of the plate was 
assumed to be uniform, because of the large ther- 
mal diffusivity Dt ~ 2 x 10 4 cm 2 /s. of the copper 
plate in the experiments. Accordingly simulations 
were carried out by the numerical method described 
in ref. .6] with a homogeneous time-dependent tem- 
perature random fluctuation of given rms amplitude 
imposed on the upper plate. This implementation 
consisted in adding or subtracting randomly temper- 
ature spikes T t at the time t with a programmed rms 
amplitude at steps separated by 0.02 s. This inter- 
val is much larger than the estimated relaxation time 
of the top plate over a distance 2L, approximately 
the wavelength of convection roll pair. Values of the 
variance A = */< (T t — < T t >) 2 > were chosen be- 
tween and 40 /iK. The range of the A values was 
taken well beyond the estimated fluctuation rms am- 
plitude during the experiments |l| of rs 1/i-ftT/v 'Hz. 
Three representative curves with 0, 3 and 40 /iK are 
shown in Fig. la) by dashed lines for e = 0.2 for q= 
2.16xl0~ 7 W/cm 2 , L = 0.106cm and T = 5.1. For 
this value of q, the calculation by El Khouri and 
Carles |l2| give t ins tab = 6.3 s and AT(t instab ) = 75 
/iK. In the simulation without imposed noise, the 
start of convection has therefore been considerably 
delayed relative to <i ns tab- The injection of random 
noise has a significant effect in developing convec- 
tion at an earlier time. In Fig. la) the three curves 
are also compared with the experimental one, shown 
by a solid line. Here we have not incorporated into 
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the simulations the delay affecting the experimental 
temperature recording, so that they could be inter- 
compared more readily, and also with predictions 
However this operation will be presented in Fig. 4. 
Further simulations with added random noise were 
carried out for e = 0.2 and 0.05 where the AT(t) 
time profiles are not shown here. 

Fig. 2a) shows a plot of the time of the developed 
convection, represented by i p , versus the random 
rms amplitude A for three series of simulations, all 
taken for a cell with T = 5.1. They are a) and b) 
e = 0.2, q = 2.16 and 3.89x 10- 7 W/cm 2 , and c) e 
= 0.05, q = 60 nW/cm 2 , ([Ra co ™ - Ra c ] = 635, 
1740 and 4200). The simulation results, shown by 
solid circles, are compared with the experimentally 
observed t p shown by horizontally dot-dashed lines. 
It can be clearly seen that noise imposition, which 
creates a vertical disturbance across the fluid layer, 
reduces the time of convection development. While 
the decrease of t p is strong for small values of A, it 
saturates at a certain level of noise amplitude. The 
gap between simulations and experiment increases 
with a decrease of [i?a corr — Ra c ) , namely as the fluid 
stability point is approached. A "critical slowing 
down" is seen in the effectiveness of the perturba- 
tions in triggering the instability. Hence this mode 
of noise introduction fails, because its amplitude is 
limited to the vertical z direction and it evidently 
couples only weakly into the convective motion. 

In parallel with the present experiments, S. 
Amiroudine 13] also carried out a systematic study 
of simulations on supercritical 3 He in a RB cell for 
several values of e and q. He used a numerical scheme 
based on the exact Navier Stokes equation as de- 
scribed in ref.Q. The resulting profiles AT(t) could 
be compared with those from experiments done un- 
der nearly the same conditions. In his simulations, 
homogeneous temperature random noise was again 
imposed on the top plate. The shift in i p showed 
less systematic trends than in the results described 
in this report. However for the same values of e and 
q as those reported above, and at zero noise, the t p 
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FIG. 1: a) the temperature profile AT(t) from exper- 
iments (solid line with noise) and from several simula- 
tions (dashed lines) at e = 0.2, q= 2.16 xl0~ 7 W/cm 2 . 
In the simulations, T—5 and uniform temperature noise 
has been imposed on the top plate with variance A(fiK) 
— 0, 3 and 40, as described in the text, b) Tempera- 
ture profile AT(t) from several simulations at e = 0.2, 
q= 2.16xl0 -7 W/cm 2 , T = 8 and imposed lateral pe- 
riodic, time independent temperature variations on the 
top plate with period 2L and amplitude B(p,K) = 0, 0.5 
and 1.5. 

values tended to be somewhat smaller than in the 
results of Fig 2a). 

Here we mention that the onset of convection in 
the simulations is further influenced by the aspect 
ratio r. The simulations described above, but with- 
out noise, were carried out in a cell T = 5.1 hav- 
ing periodic lateral boundaries. Further simulations 
with zero noise for e = 0.2 with T = 8.0, 10.2, 20.5 
and 41.0 were carried out, and showed a decrease of 
the convection development time from w 90 s, tend- 
ing to a constant value of w 60 s. above T = 20. 
This shift in the onset of instability is due to the 
decreased finite size effect which the rising plumes 
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FIG. 2: a) The time for effective development of convec- 
tion, characterized by t p , versus A (homogeneous tem- 
perature noise imposed on the top plate). The horizon- 
tal dot-dashed lines indicate the observed t p , corrected 
for instrumental recording delay, b) The time for effec- 
tive development of the convection, labeled by t p versus 
B (lateral time-independent periodic temperature vari- 
ations). The horizontal dashed line indicates t p as ob- 
tained by experiment, corrected as before. 

experience with increasing T, in spite of the periodic 
boundary conditions. This can be seen by compar- 
ing the curves labeled "O" in Figs la and lb with T 
= 5 and 8 respectively. 

The next step in our attempts, stimulated by 
communications with P. Carles, was introducing 
perturbations into the simulations via some time- 
independent lateral variation proportional to sin 
(2ttx/P) where P is the period. We opted to intro- 
duce again a temperature variation in the top plate 
with an amplitude B (in /xK) and period P=2L, 
nearly the same as the wavelength of a pair of con- 
vection rolls. The temperature of the bottom plate 
was kept homogeneous. This "Gedanken Experi- 
ment" implies that the material of the top plate 
permitted a temperature inhomogeneity, which of 
course is not realized in the experiment. However 
a small lateral temperature excursion can trigger 
the same kind of non-homogeneous perturbations as 
those which, in the real experiment, provoke the on- 
set of convection. One possible origin of such pertur- 



bations, besides thermal noise, could be the rough- 
ness of the plates or their slight deviation from paral- 
lelism. Such geometrical defects could of course not 
be implemented in the numerical simulations with 
the meshsize used, which is why we elected to force 
a small temperature perturbation instead, with sim- 
ilar effects on the onset. As a control experiment, 
we also made a simulation with P = L. 

Fig lb) shows representative profiles AT(t) for the 
parameters e = 0.2 and q= 2.16xl0~ 7 W/cm 2 and 
with B = 0, 0.5 and 1.5 fiK, and for r=8. As B 
is increased from zero, there is a large decrease in 
the time for convection development, represented by 
t P) which is plotted versus B in Fig 2b). The hor- 
izontal dashed line shows the t p from the experi- 
ment, and this plot is to be compared with Fig. 2a). 
For an inhomogeneity amplitude of only B= 0.5/iK, 
t p is nearly the same for simulations and experi- 
ment. By contrast, simulations with B=2/iK and 
P=L (not presented here) show no difference from 
those with B=0. Hence the nucleation of the con- 
vection is accelerated if the period is in approximate 
resonance with the wavelength of a convection roll 
pair. The values of steady-state AT and t osc are 
only marginally affected by the noise. 

We note from Fig. lb) that the simulation curve 
calculated for B = shows the fluid not convert- 
ing until w 70 s. For the curves with B= 0.5 fiK., 
the start of deviations from the stable fluid curve 
cannot be estimated well from Fig. lb) but is readily 
obtained from the data files, which tabulate AT(i) 
to within 1 nK. For B = 0.5 fiK, systematic devia- 
tions SAT(t, B) = [AT(t,B=0) - AT(*,B)] increase 
rapidly from 1 nK for t >8 s (where AT «85 //K), 
a value comparable with the predicted imstat>= 6.3 
s., AT(ti nsta b) = 75^iK^^. However a comparison 
with predictions becomes more uncertain as B is in- 
creased and no longer negligible compared with the 
steady-state AT. Then it is expected that the base 
Piston-Effect heat flow will become itself influenced 
by the perturbations. In that case the stability anal- 
ysis in 0| becomes irrelevant, since the base flow, 
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FIG. 3: Snapshots in 2D for the RB cell with an aspect 
ratio of 8, of simulations with B = 0.5 fiK at various 
times t after starting the heat current q. The temper- 
ature contour lines and their evolution are described in 
the text. At the time t = 9.9s, the fluid instability has 
just started near the top. 

the stability of which is analysed, has been signifi- 
cantly modified by the perturbations. We also note 
that the time interval St = [t p — i; ns tab] between the 
first sign of instability (SAT > 0) and t p is » 20 s, 
and roughly independent of B. This represents ap- 
proximately the period taken by the convection to 
develop and for the plumes to reach the top plate 
boundary. 

In Fig. 3, we present a series of 2D "snapshots" at 
various times for the simulation with B = 0.5/xK, 
showing the temperature contour lines (in color) for 
the RB cell. The "warm" side is shown by red, 
T(t, z = 0) and the "cold" side by mauve, T(z = L) 
= const. At t= 8 s. the fluid instability has just 
started near the top of the layer, while near the bot- 
tom the temperature contour lines are still horizon- 
tal. At t = 27 s., where the peak of AT(t) at z = L 
has been reached, the warm plumes have reached 
the top plate, and the "cold" piston effect is about 
to start, causing the bulk fluid temperature to drop 
and AT(t) to decrease. The transient process con- 
tinues with damped oscillations of AT(t). Steady 
state convection is reached at t= 80s, with a pair 
of convection rolls having a wavelength of « 2L, as 
expected. 

In Fig. 4 we show the profiles AT(t) from the ex- 
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FIG. 4: Comparison of the profile AT(t) from experi- 
ment and from simulations with B = 0.5 \iK . To make 
the comparison realistic, the simulations have been con- 
voluted with the same "instrumental" delay time r —1.3 
s. which has influenced the shape of the experimental 
curve. 

periment and from the simulations with a periodic 
perturbation amplitude B = 0.5 /iK. For an opti- 
mal comparison, the delay affecting the experimen- 
tal temperature recording was incorporated into the 
simulation curve. For this, the delay function with 
the instrumental time constant r = 1.3 s. |]| was 
folded into the simulation curve by a convolution 
method. This operation retards the initial rise of 
the temperature drop by the order of 2-3 seconds, 
and brings both experiment and simulations into fair 
agreement in the regime where the fluid is stable. 
The time t p for the maximum is now closely the 
same for both experiments and simulations. How- 
ever beyond the predicted instability time tinstab = 
6.3 s., the experimental curve starts to deviate more 
rapidly with time than do the numerical simulations 
from the calculated curve for the fluid in the sta- 
ble regime. As discussed previously Q, for these pa- 
rameters of e and q the experiment does not show 
damped oscillations, which are observed for higher 
values of q. In the steady-state, the agreement is 
very good. 

Our goal has been to show that injecting a small 
temperature perturbation into the top plate, pro- 
duces for the simulations an earlier start in the con- 
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vective instability, which becomes consistent with 
experimental observations. For this, we have lim- 
ited ourselves to an example at a low value of 
[Ra c °" — Ra c ], where the delay has been particu- 
larly large with respect to the experiment. 

SUMMARY AND CONCLUSION 

We have presented a comparison of numerical sim- 
ulations with experimental data investigating the 
transient to steady convection after the start of a 
heat current through a supercritical 3 He layer in a 
RB cell. Here the temperature drop AT(t) across 
the fluid layer versus time t was studied. The aim 
was to understand and to reduce the discrepancy 
between experiment and simulations in the time of 
the convection development, as detected by AT(t). 
Simulations for one set of fluid parameters (where 
the largest discrepancy had been observed) are re- 
ported with imposed temperature variations on the 
top plate. Satisfactory results were obtained for spa- 
tial lateral temperature variations with an amplitude 
of 0.5 /iK and a period approximately equal to that 
of the wavelength of a convection roll pair. As the 
perturbation amplitude is further increased, the de- 
velopment of convection occurs earlier than the ob- 
served one. 
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